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Abstract. We investigate stationary solutions of a thin-film model for liquid two-layer flows in 
an energetic formulation that is motivated by its gradient flow structure. The goal is to achieve a 
rigorous understanding of the contact-angle conditions for such two-layer systems. We pursue this 
by investigating a corresponding energy that favors the upper liquid to dewet from the lower liquid 
substrate, leaving behind a layer of thickness h t . After proving existence of stationary solutions for 
the resulting system of thin-film equations we focus on the limit h* — > via matched asymptotic 
analysis. This yields a corresponding sharp-interface model and a matched asymptotic solution that 
includes logarithmic switch-back terms. We compare this with results obtained using T- convergence, 
where we establish existence and uniqueness of energetic minimizcrs in that limit. 

1. Introduction. Understanding stability and dewetting behaviour of thin liq- 
uid films coating a solid or a liquid substrate is important in many technological 
applications and natural phenomena on the micro- to nano scale. They range from 
tear films of the human eye to organic photovoltaics to numerous applications in the 
polymer based semiconductor industry. 

Typical film thicknesses for these applications may range from tens to hundreds 
of nanometers and, depending on the material composition, may be susceptable to 
rupture and formation of holes due to intermolecular forces. Such rupture processes 
typically initiate complex dewetting scenarios, where holes grow further and their 
trailing rims merge into polygonal networks which eventually decay into patterns of 
droplets, that evolve on a slow time scale towards a global minimal energy state. 

The present study focusses on liquid substrates, that energetically favor an inter- 
face with the underlying solid. In this case the stages of the dewetting process for the 
upper liquid proceed to some extend in parallel to those exhibited during dewetting 
of a liquid film from a solid substrate. The latter system has been investigated much 
more intensely in recent decades, both, experimentally and theoretically. Examples 
of the complex pattern formation can be found in Sharma et al. [l[ or Seemann et al. 
and further experimental and theoretical investigations in numerous references in 
the recent reviews by Craster & Matar Q and Herminghaus et al. jij. 

For liquid-liquid dewetting, experimental studies depicting some of these dewet- 
ting stages have been conducted by various groups such as by Scgalman & Green [f| , 
Lambooy et al. @, Slep et al. Q or Wang et al. [|| for the standard system of liquid 
polystyrene (PS) on a liquid polymethylmethacrylate(PMMA) substrate. They in- 
clude investigations of rupture and hole growth, dewetting dynamics and equilibrium 
contact angles the liquid droplets make with the underlying liquid layer, where now 
the contact line is fixed by two angle conditions instead of one, i.e. Youngs law is 
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replaced by the Neumann triangle construction |9|. Following the pioneering study 
by Brochard-Wyart et al. [l(| , where various dewetting regimes are derived and anal- 
ysed, stability of liquid-liquid systems were investigated by Danov et al. , Pototsky 



et al. [12|, Golovin & Fisher [13(. Stationary states and the dynamics towards sta- 



tionary states were studied by Pototsky et al. 14|, by Craster & Matar and by 



Bandyopadhyay & Sharma [iff and for the case of gravity-driven liquid droplets on 
a inclined liquid substrate by Kriegsmann & Miksis hj. 

Interestingly, direct quantitative comparisons of theoretical with experimental 
results, in particular on micro- and nano-scale, regarding for example the morphology 
of the interfaces, as performed in Kostourou et al. [IH, or equilibrium values of the 
Neumann triangle, as discussed in Q , still leave many issues in need to be explained, 
such as the dependency of the morphology of the interfaces on the rheology of the 
liquids, their layer thicknesses or material parameters. On the other side, even for the 
simplest mathematical models of Newtonian two-layer liquid systems, mathematical 
theory is still largely open and this is the main focus of the present study. 

Here, we are guided by the many similarities to dewetting from a solid substrate, 
and expect that some of the mathematical analysis developed for liquid films dewetting 
from a solid substrate can be carried over to liquid substrates. In particular, we aim 
to extend the existing theory for liquid droplets on a solid substrate to the situation 
on a liquid substrate. 

As a starting point we recall the work by Bertozzi et al. [l9[ regarding stationary 
states and coarsening of droplets on solid substrates, where they showed existence of 
smooth global solutions for positive data with bounded energy for the no-slip lubrica- 
tion equation. In addition they prove existence of global minimizers and determine a 
family of positive periodic solutions for admissible intermolecular potentials consisting 
of long-range attractive and short range repulsive contribution, and investigate their 
linear stability. Further extensions were given in Laugesen & Pugh (20| . where linear 
stability of stationary solutions for the thin film equation with Neumann boundary 
conditions or periodic boundary conditions was investigated. Extensions of the ex- 
istence theory to thin film equations accounting for slip at the liquid-solid interface 
were given in Rump et al. [2l[ and Kitavtsev et al. [22[ • Convergence to stationary so- 
lutions of the one-dimensional thin-film equation and the number of stationary states 
was recently investigated by Zhang [23 1 . 

The extension of the existence theory to two-layer liquid systems is given in Sec- 
tion 3, after the formulation of the problem. With the appropriate energy functional 
for the two-layer system of coupled thin-film equations for the interfacial heights hi 
and h 2 (see sketch in figure [TTTj) together with Neumann boundary conditions, we show 
existence of smooth stationary solutions as well as existence of a global minimizer for 
the steady state problem. 

Admissible intermolecular potentials for the liquid-liquid system are of the form 



<fi(h 2 - hi) 



h 2 — hi 



h 2 - hi 



(1.1) 



where hi is the height of the liquid-liquid interface, h 2 the height of the free surface 
and its minimal value < is attained at h t . We note that such potentials are 
widely used in the literature, see e.g. the review 24|. A choice that is related to the 
standard Lennard- Jones potential and typically used in experiments is (n, £) — (2, 8), 
see e.g. @- 
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Fig. 1.1. Sketch of liquid droplet with surface tension 0-2 between air and liquids on top of a 
liquid layer with interfacial tension o\ between liquid^ and liquid i. 

The Neumann triangle construction for contact angles at a triple junction is 
thereby replaced by properties of approximate contact angles resulting from the par- 
ticular structure of the surface free energy <j). Starting from this energetic formulation 
we seek to understand in Section 4 and 5 how the Neumann triangle construction is 
attained as a limit h* — > 0. This in mind, we first derive in Section 4 a sharp-interface 
model in the limit of h* — > using matched asymptotic analysis. This yields the 
appropriate expression for the equilibrium contact angles of these droplet solutions 
and as a result the corresponding Neumann triangle. 

We find that, while the equilibrium contact angle is easy to obtain, as expected, 
the complete asymptotic argument needs to include logarithmic switch-back terms. 
We note, in retrospect, since equilibrium droplet solutions for solid substrates can be 
considered as limiting cases for liquid lenses, similar terms should also appear in the 
matched asymptotic derivation for that problem. 

Finally, in Section 5 we rigorously show existence and uniqueness of the limit 
/i* — > within the framework of T-convcrgcnce. We obtain a sharp-interface problem 
for which we compute the Euler-Lagrange equations, from which one can immediately 
read off the contact angles. Existence and uniqueness of minimizers is shown here 
using a rearrangement inequality 

2. Formulation. We consider a layered system of two immiscible Newtonian 
liquids with negative spreading coefficient 0*. We assume the layered system lives in 
the two phases fii and defined by 

fi x (f) = {(x,y,z) £l 3 :0<K hi(t,x,y)}, (2.1a) 
iJ a (t) = {(*, V, z)GR 3 : ht(t, x,y)<z< h 2 (t, x, y)}, (2.1b) 

for all t > as sketched (in 2D) in figure H~T1 Typical applications use liquids such 
as PMMA for the liquid substrate fii and PS as the upper liquid fl2, both on a 
scale where gravity can be neglected and, for this study, unentangled and density 
matched. These simplifying assumptions allow us that the flow of the viscous and 
incompressible liquids in each phase fli (i = 1, 2) is governed by the Stokes equation 
and the continuity equation: 



-V • (-p 4 I + Mi(Vui + VuJ)) = f; , 
V • Ui = 0, 



(2.2a) 
(2.2b) 



4 



Jachalski et al. 



together with a kinematic condition at each free boundaries z = h i7 i.e. 



{e z d t h t - Uj) • n.; = 0. 



(2.3) 



Here, n; denotes the outer normal. 

At the solid substrate a no-slip and an impermeability condition are imposed. At 
the liquid-liquid interface the velocity is continuous, i.e. u 2 = Ui, and the tangential 
stress is continuous and the normal stress makes a jump proportional to the mean 
curvature [[— pjI + /ij(Vuj+Vu^")]]i i2 ni = 2(T li2 Kn 1 . The dewetting process is driven 
by the intermolecular potential of the upper liquid layer, i.e. f 2 = — V0. We assume 
that the thickness of the lower layer is sufficiently thick, so other contributions to the 
intermolecular potential are negligible, i.e. fi = 0. 

In addition we assume that the the ratio eg = H/L of the vertical to horizontal 
length scales is always small and the twodayer system can be approximated by a 
thin-film model. We denote by L the length scale of the typical horizontal width and 
H = /imax the maximum of the difference h% — h\. 

Detailed derivations of thin-film models for liquid-liquid systems are given for 
example in (l7j . fl2j or recently, by accountung for interfacial slip, in (2a |. For 
the convenience of the reader we note here the choice of non-dimensional variable 
and parameters, where the non-dimensional horizontal and vertical coordinates are 
given by x = x/L, y = y/L and z = z/h max , respectively, and h = h/h max . The 
non-dimensional pressures pi = Pi/P and the derivative of the non-dimensional in- 
termolecular potential (f>' e = 4> jP are scaled such that P = <fr*n£/((£ — n)/i max ) and 



attains the minimal value min0 e = — 1 at h = e, where e = h*/h max is the non- 
dimensional thickness of the ultra-thin film. For the dynamic problem, the velocities 
are scaled with the characteristic horizontal velocity of the dewetting upper layer, 
such that for the non-dimensional horizontal and vertical velocities we have u = u/U, 
v = v/U and w = w/egU, respectively, with U = £^112/1^2, and the non-dimensional 
time t= (U/L)t. 

For the remainder of this paper it is convenient to introduce the ratios a = o\ j 02 
and /i = /X1//X2 of surface tensions and viscosities, respectively, and drop all the 
"~" . Within this approximation the normal and tangential stress conditions at the 
free surface /i 2 and at the free liquid-liquid interface h\ yield the expressions for the 
pressures p 2 and pi 



respectively. Under these assumptions the coupled system of nonlinear fourth order 
partial differential equations for the profiles of the free surfaces hi and /i 2 takes the 
form 



hence 




(2.4) 



pi = -aAhi - 4>' s {h2 - hi), 



p 2 = -A/i 2 + <t>' e {h 2 - hi), 



(2.5) 



d t h = V • (Q ■ V P ) , 



(2.6) 



where h = (hi , h 2 ) T is the vector of liquid-liquid interface profile and liquid-air surface 
profile. The components of the vector p = (pi,p 2 ) T are the interfacial pressures given 
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in (|2.5[) . The gradient of the pressure vector is multiplied by the mobility matrix Q 
which is given by 



Q = - 



M 

3 

hi hl{h 2 -h x ) 
3 2 



h\ hljhi-hx) 
3 2 

(h 2 - h x ) 3 + fti/i 2 (/i 2 - fti) + 



3 J 



(2.7) 



We proceed to first show existence of stationary solution to the system (|2.4[) - (|2.7p 
with Neumann and Dirichlet boundary conditions on a finite domain = (0, L) Cl. 

3. Energy functionals and existence of stationary solutions. Consider 
the two- layer thin film equations (|2.5[) - (|2.7[) defined on = (0, L) c K with 
Neumann boundary conditions: 



d x hi = d x h 2 = d xxx hi = d xxx h 2 = for x £ dfl. 



(3.1) 



The energy functional associated to the gradient flow of the lubrication equation is 
given by 



EeihxM) 



^\d x hx\ 2 + ^\d x h 2 \ 2 + ^2-^) 



dx 



(3.2) 



where the potential function <j) e is given as in (|2.4p with (n,£) = (2,8). The relation 
to the thin-film equations is Pi = 8E e /8hi. From (|2.5[) - (|2.7[) and (|3.1|) conservation of 
mass follows 



(3.3a) 
(3.3b) 



hx(t, x) dx = mi, 
{h 2 {t,x) — hi(t,x)) dx = m 2 for all t > 0, 



where mi and m 2 are positive constants. Any stationary solution of (|2.5[) - ([2.7p with 
(|3.1[) satisfies 



= Q • d x p 



(3.4) 



in J7, where the mobility matrix Q is not singular i.e. dctQ 7^ for all hi, hi — h\ > 0. 
Therefore, one obtains that any positive stationary solution of (|2.5[) - (|2.7[) satisfies 
dxPi = d x p2 = in Q. This in turn is equivalent to 



(rd xx hi = ~(j)' € {h 2 - hi) - A 2 + Ai, 
d xx h 2 = (/)' e (h 2 - hi) - Ai, 



(3.5a) 
(3.5b) 



where constants A2 and Ai are Lagrange multipliers associated with conservation of 
mass (|3.3a[) and (|3.3b|) . respectively. To solve (|3.5p let us consider the equation for 
the difference 



h(t,x) = hi{t,x) — hi(t,x) 



which reads as follows 



d xx h = G + (f>' e (h) + -A 2 
a a 



(J + 1 



(3.6) 
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For brevity we set 



1 CT+1 

-A 2 H Ai. 



According to |19| for positive P, there exists a so called droplet solution h to (|3.6p 
satisfying boundary conditions (|3.ip . such that fay + L/2) is an even function and 
monotone decreasing for y <E (0,L/2). 

For this solution the asymptotics and main properties are derived in the next 
section, here we consider h as a known analytical function and integrate (|3.5p two 
times w.r.t. x. We then obtain a solution to (|3.5[) with (|3.ip in the form 

fci = \- h - l^-x 2 + Cx + d, 

(7 + 1 2 a + 1 

fa = ~^-rh- \^-x 2 + Cx + Ci. (3.7) 
(7 + 1 2 a + 1 

Using now again (|3.1|) one obtains that A2 = 0, C — and 

/11 = -ft + Ci, 

CT+ 1 

fa = -^—h + C x . (3.8) 

cr + 1 

The additive constant C\ and the remaining Lagrange multiplier are determined from 
the conservation of masses p.3ap and (|3.3b[) , respectively. We conclude that the so- 
lution (|3.8p is given by combination of two symmetric droplets with constant outer 
layer. The next theorem establishes existence of a global minimizer to the energy 
functional (13.21) and shows that it satisfies (13.51) with (13.11). 



Theorem 3.1. Let il be a bounded domain of class C ' 1 in M. d , d > 1 and let 
m = {mi, mi) with mi,mi > 0. Then a global minimizer of E e (-,-) defined in i3.!2\) 
exists in the class 

X m := \{h x M) G H\nf :mi = f h u m 2 = [fa -fa), h 2 > fa\ , (3.9) 
I Jn Jn J 

For d = 1 and = (0, L) the function fa — hi is strictly positive and (fa, fa) are 
smooth solutions to the ODE system VS. 5\) with H3.1\) and 

Ai = - / (f>' e {fa - fa) dx, A 2 = 0. (3.10) 

Li .In 



Proof. Even though the proof proceeds very analogously to the one of Theorem 
2-3 in [19( using direct methods of the calculus of variations, for the convenience of 
the reader we give here for our system. 

Let (hi, hfykeN be a minimizing sequence which exists since E e (mi, m 2 + mi) < 
00. Observe that 4> t (-) is bounded from below by a constant. Hence, a constant Ci 
exists such that 



/ \d x h\ | 2 + \d x h\^ dx < C 2 for all k e N. 
Jn 



(3.11) 
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Rcllich's compactness theorem implies that a subsequence (again denoted by (h\, /i^fceN) 
exists which converges strongly in L 2 (f2) 2 and pointwise almost everywhere to (hi, h 2 ) £ 
iJ 1 (il) 2 . By Fatou's lemma, we deduce that also <f> e (h 2 — hi) lies in Using the 

weak lower scmicontinuity of the norm, we obtain 

E e (hi,h 2 ) < lim inf E e (h k ,h%) = inf E e (h 1 ,h 2 ). 

fc->oo (h 1 ,h 2 )£V 

Consequently, (hi, h 2 ) is a minimizcr and the intcgrability of 4> e (h 2 — hi) implies that 
h 2 — hi > almost everywhere in Q. 

Let now d = 1 and f2 = (0, L). The estimate p. lip implies 



\d x (h k ~h k )\ 2 dx<C 3 . 



u x[ r k ^' 2 

'O 

Next, the boundedness of J 4> e (h k ~ h k ) dx and definition of 4> t imply a uniform in k 
bound on \\(h k ~ Tt.^ ) 4 1 1 2 - Using this, one can estimate 

<3||a»(/>;-/ ! ;)|| 2 ||(/i|-/!f)- 4 ||2<c 1 , 

where C4 is constant. Owing to the continuous embedding of H /1,3 (0, 1) into L°°(0, 1) 
the strong positivity of h 2 — hi follows. This in turn implies the differentiability of the 
function F t (s) := E e (h\ + sipi,h 2 + sip 2 ) considered with fixed (ifi,ip 2 ) e i? 1 (0,L) 2 
for sufficiently small s. Since (hi,h 2 ) is a minimizcr, by differentiation of F £ (s) at 
s = 0, we obtain that 

L 

(-vd xx hi - 4>' e {h 2 - hi))cpi + (-d xx h 2 + (f>' e (h 2 - hi))if 2 dx = 0, 
for all ((pi,(p 2 ) e i/ 1 (0,i) 2 such that 

L pL 

ipi dx = / (if 2 ~ <Pi) dx = 0. 
Jo 

Without this constraint, using Lagrangian multipliers and testing with general (-01 , i/) 2 ) G 
ff^O,^) 2 this yields 



(-<rd xx hi - d xx h 2 )i\)i + (-d xx h 2 + 4>' £ )ip2 dx — — I I <f> e dyip 2 dx = 0. 



L JQ JQ 



Standard elliptic regularity theory then implies that (hi, h 2 ) are smooth solutions to 
(j3~5j) together with {(37TJ) and (pTTT>|) . □ 



Remark 3.2. Note, that for Dirichlet boundary conditions we can proceed as 
follows: Let us impose on system i3. 5\) the Dirichlet boundary conditions 



hi(0) =hi(L) = A, h 2 (0) = h 2 (L) = B, (3.12) 

such that 

B-A= min h(x), (3.13) 
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where h is the Neumann solution to \3. 6}) defined above. In this case it follows 
again that h 2 — hi = h. Consequently, hi and h 2 are given by \3. 7| ) with constants 
Xi,X 2 ,C, Ci determined uniquely by liS.lty) and conservation of mass A3. 3a\ )- l3.3b\) . 
Using i3.13\) and the asymptotics for h one obtains that the leading order of the solu- 
tion \3. 7| ) as e — > has now the form 

h 2 (x) = -±\ 1 ((x-§) 2 -s 2 )+C 1 , 
hi{x) = h 2 (x) = -1-^ ((x - |) 2 - s 2 ) + Cx, 
where uj = (L/2 — s, L/2 + s) and 

(A 2 - (a + l)Xi) 2 ' 

In contrast to the solution of (|3.14j) with Neumann conditions, solutions with Dirichlet 
boundary conditions are not constant but quadratic in the ultra-thin layer (0,L) \cj. 

4. Matched asymptotic solution and contact angles. Note first that the 
system of equations for hi and h 2 (|3.4p is equivalent to following system for hi and h 



x e uj 
x e uj 

x€{0,L)\w, (3.14) 



= d x (-ad xx hi - 4>' £ {h)) , 
= 3 X ( ^—d xx h + cp' £ (h) 

\ (7+1 



(4.1a) 
(4.1b) 



where we denote a = ai/cr 2 , see [25[ for a more detailed derivation. For our asymptotic 
analysis, this is convenient, since now for the variable h = h 2 — hi we can distinguish 
the core droplet region, which we will call the "outer region" and the adjacent thin 
regions of thickness e, which we call the "inner region". We will derive a sharp- 
interface limit using matched asymptotic analysis in the limit as e — > 0. For this 
we first write the equations in the form such that the intcrmolecular potential is 
small in the core region and becomes order one in the adjacent thin regions. Using 
(n,£) = (2,8) we define 



(4.2) 



4.1. Stationary solution for h. As we will later show rigorously, we can as- 
sume that the droplet is (axi) symmetric and the profile a decreasing function, so that 
without loss of generality, the maximum of h is at the origin of our coordinate system. 
Now consider the problem h and x > 



= d x 



-^-d xx h-U> ( h - 

a + 1 £ \ £ 



(4.3a) 



lim h = /loo, hm d x h = 0, lim d xx h = 0. (4.3b) 

X— >OG X— >00 X— fOO 
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We can integrate this twice and use the conditions as x — > oo to fix the integration 
constants to obtain 

W - ^±I^(i)-.(^)-I V (^) (4.4, 

The solution to this problem shows a steep decline in height towards 0(e) in an e-strip 
around x = s, where we would like to determine the apparent contact-angle. This can 
be obtained by writing the problem in so-called outer and inner coordinates, valid in 
the core and the adjacent thin regions, and matching as e — > 0. Interestingly, while 
it is easy to obtain the condition for the contact angle, it turns out that in order to 
carry out the complete matching consistently, we need to go up to second order in the 
matching, in order to account for the logarithmic switch-back terms, that come into 
play in this problem, see j2(| for a discussion of these terms. 

Note, that the coefficient (<r + l)/er can be removed by rescaling x appropriately, 
leading to the classical problem of a droplet of height h on a solid substrate. Interest- 
ingly, to our knowledge for this problem the above mentioned logarithmic switch-back 
terms have not been noticed before. 

Outer problem. The symmetry of the problem leads us to the condition that 
at the symmetry axis x = we have 

d x h = (4.5) 

It is also convenient to normalize the height such that h(0) = 1. In this case we obtain 
from (|4.4p an algebraic equation for hoo and e that can be approximated as e — > 0. 

„_.(!) (£)(!-*_) (4.0, 

Solving this by making the ansatz for the asymptotic expansion for 

hoc = eh xfi + e 2 h ooA + £ 3 h 00)2 + 0{e i ) (4.7) 

we obtain 

1 45 
^oo,o = lj hooi = — , /ioo. 2 = 777:7 (4-8) 
lb 512 

Next, we assume that the asymptotic solution to the outer problem can be represented 
by the expansion 

h(x;e) = fo{x)+eMx)+e 2 f 2 (x) + 0(e 3 ). (4.9) 
The leading order outer problem then becomes 



which has the solution 



d x fo = -^±1(1-/0) 
/o(0) = 1 

/o(,) = ~^ 2 + l. 
16 a 



(4.10a) 
(4.10b) 

(4.11) 
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Hence, the leading order outer solution will vanish as x approaches the location 

However, the full solution docs not vanish and will be obtained by matching to the 
solution of the "inner" problem near s. We will find that in order to complete the 
solution, we will need to solve the expansion up-to second order. We find for j\ and 

h 

d x h= f -±-^°-±±x, (4.13a) 
x 16 cr 

A(0) = (4.13b) 

and 

d x h = ^- + - f ! Jj + i ( /o - 1) - | ( /o + 3 ) + 2h) + l—x (4.14a) 

x x \6 Jo Id 3 J o a 

MO) = (4.14b) 

Inner problem. The solution of the inner problem lives in an e neighborhood 
of x = s and extends towards x — > +oo. It will be matched to the outer problem 
in the other direction. Hence we introduce the inner variables v(z) and independent 
variable z via 

h(x) = ev(z; e) and x = s + ez. (4-15) 
Rewriting the problem (|4.4j) in these coordinates and making the ansatz 

v(z;e) = v (z) + sv 1 (z) + O(e 2 ) (4.16) 
we find to leading order the problem 



a+l / 1 . 1 

2vk 



d z v = - x j2— x \-^ + -g + - (4.17) 



and to O(e) the problem 



civ - J? a + 1 I 1 4- 1 + 3 3^gK-l)+8^(l-^) 

d zVl - -^2—^-^ + ^ + - 2wo(4u 6_ 3w 8_ 1) 



We solve and match in the inner coordinates and obtain by expanding v at z = — oo 



y/3 1(7 + 1 4 I a 1 

vq + ev\ = — — \ z + a 1 - -4/3- 



2V cr 9\j a + l z 

( 3 a + l 2 73(7 + 1, 
- £ -__ z =' + __( ai -l) z -ln(-z) 



+c nnV^TT (ai + 1) I + -»- 
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For the corresponding outer expansion we have 



fo + efi + e 2 h 



\/3 fa+l 



z 



+£ 



a 

3 cr + 1 
16 a 



4 / a 1 

1 - -,/3 

9V (J+l z 



z — ln(— z) 



19 

- B -M.) 



(4.20) 



We note that all of the terms in the first row of (|4.19[) and (|4.20[) match provided 
a\ = — 1. The terms in the second row match provided 



C 



(4.21) 



where we note the appearance of a so-called "logarithmic switch-back" term ln(e). 
Hence, the composite solution is 



i'O 



\/3 j a + 1 x — 
~2 



(4.22) 



16 cr + 1 



V3CT + 1 
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\/3 /cr + 1 



with s given in equation (|4.12[) and for x < s. For x > s only the inner expansion 
h = evq + e 2 v\ remains. 

Stationary solution for hi and h<x. To complete the solution, we determine the 
solution to the liquid- liquid interface h\ simply by adding equations (|4.1ap and (|4. lb|) . 
Then we integrate thrice and use the far-field conditions d xx h,d x h,d xx h\,d x h\ — > 0, 
h — > hao and hi — > d as x — > ±oo to fix the constants. This results in hi being 



h i 



1 



(7+1 



and equivalently hi is 



hi 



<J + 1 



(h — hoo) + d + h c 



(4.23) 



(4.24) 



Here d denotes the height hi as x — » ±oo and we assume d to be large enough so 
that hi never becomes negative. Also note that for other boundary conditions, such 
as e.g. Dirichlet conditions mentioned in the previous section, further contributions 
will arise. Unlike the solutions for droplets on solid substrates, here new families 
of profiles for the ultra-thin film connecting a droplet to the boundaries or to other 
droplets arise. 
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Fig. 4.1. Top: Comparison of the composite solution h c (dashed curve) with the numerical 
solution of l|4,3jl (solid curve), for e = 0.2 and a = 1.2. Bottom: Solutions for hi and hi recon- 
structed from the composite solution h for a = 1.2 and e = 0.05 (dashed curve) and e = 0.2 (dotted 
curve) and the solution to the sharp interface model l|4.25|l (solid curve). 



In figure 14.11 we compare our asymptotic solution with the numerical solution. 
Observe that the 0(e) solution already gives an excellent approximation of the exact 
solution for e = 0.2, where the exact solution is approximated by the higher-order nu- 
merical solution of the boundary value problem. This suggests that a sharp-interface 
model should also be a good approximation of the full model. The sharp-interface 
model for h is simply the leading order outer problem for h, but now with boundary 
conditions, that result from the leading order matching. Hence, we obtain 



— dxxxfo 



(4.25a) 



with boundary conditions 



/o = and d x f = d z v 



-2^il$(l) at x = ±s 



(4.25b) 



where the contact angle has been determined via matching. Equivalently, for the 
half-droplet, by imposing symmetry we have the boundary conditions 



/o(s)=0, 9 K /o(s) = -y-2^^$(l) and d x f (0) = (4.26) 

In the following section we show how the sharp-interface model is obtained via 
T-convergence and also proof existence and uniqueness of its solutions. 

5. Sharp-interface limit via T-convergence. In this section we investigate 
properties of stationary solutions of the sharp-interface two-layer model. Such model 
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can be obtained by considering the limiting problem e — > in the framework of 
T-convergence. For one-layer systems corresponding minimizcrs are known as meso- 
scopic droplets [2l|. In this approach equilibrium contact angles can be directly de- 
duced from the Euler-Lagrangc equation of the resulting T-limit energy. On the other 
hand showing an equi-coerciveness property we have that the sequence of minimizcrs 
of E £ converges to a minimizer of the T-limit energy Eoo . 

For boundary conditions fa = h 2 and certain domains we show that solutions of 
the minimization problem min^ j i2 ) exist and are unique up-to translations. The 



main technique used here is the symmetric- rearrangement, see e.g. [27 , 

For the section to come we consider energies such as in (|3.2[) . For later convenience 
we define W e (h) = (4> e (h) - $(1))/|$(1)| where as before W(h) = W e (h/e) is inde- 
pendent of e. The shift by |$(1)| has the advantage of working with a non- negative 
energy without changing the Euler-Lagrange equations. 

With these definitions consider the following family of minimization problems. 
For Q C R" bounded with Lipschitz boundary and mi,iii2 > given and e > we 
look for minimizcrs of E e : X m R°° defined as 

Ee(hi, h a ) = J ||V^| 2 + ^IV^I 2 + |$(1)| W 6 (h 2 - fa) (5.1) 

with a > and W £ (h) = W(h/e) as e — > 0. Note that nonnegativity h\ > is not 
enforced because otherwise extra terms in E e are required. 

5.1. T-convergence. For a given domain il C 1™ bounded with Lipschitz 
boundary define the space of admissible interfaces as before in (|3 .9[> by 

X ia :=\(h 1 ,h 2 )€H 1 (n) 2 :mi= f h u m 2 = [fa -fa), h 2 >fa\, (5.2) 
I Jn Jn J 

and in general let W : M — > R be a function with the properties 

(i) W > and W(h) = ^ h = 1. 

(ii) W(h) ^ +oc > 1 and W(h) < 1 for h > 1. 

We want to investigate the family of minimization problems (|5.1[) . First we note 
that the sequence E e is equi-coercive in the weak topology of H 1 ^}) 2 . This is a simple 
consequence of 

E £ {fa,h 2 )>c{\\Vfa\\ 2 Hl + \\Vh 2 \\ 2 H1 ), 

which holds for all (fa , h 2 ) € X m . Together with the T-convergence the equi-coercivity 
implies the following abstract convergence result. We know that any sequence {fa^ n , h 2 ^ n } 
of minimizcrs to the energies E Sn has a weakly converging subsequence (fa, n , h 2 ,n) — 1 
(h\, h 2 ). Furthermore the limit (fa, h 2 ) is a minimizer of the T-limit energy E^. This 
relates minimizers of E s to minimizers of the T-limit E^ . 

Now we investigate the r-limit of (|5.1|) in the topology of weak convergence in 
the space _ff 1 (fi) 2 . We recall the definition of T-convergence , see also 29|, 3pj . 



Definition 5.1. We say that a sequence E £ : X —> R°° Y -converges in X in the 
weak topology to E^ : X — > R°° if for all x G X we have 
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(i) (lim-inf inequality) For every sequence {x e } C X weakly converging to x there 
holds 

Eoo(x) < liminf E e (x e ). 

e 

(ii) (lim-sup inequality) There exists a sequence {x E } C X weakly converging to 
x such that 

Eoo{x) > limsup£ , e (x e ). 

6 

The function Eoq is called the T -limit of {E e ), and we write E^ =r-lim e E e . 

The key proposition to compute the overall T-limit is to consider the T-limit of 
the potential separately. Here we use that weak convergence in H 1 implies strong 
convergence in L 2 and the right continuity of q i— > J x{h > q} for any given h g H . 

PROPOSITION 5.2. Consider the functional F E : fl' 1 (Q) R°° defined as 

' J a W e (h) h€X m , 



Fe{h) 



othe 



whe 



Then 



h = rr 



X m = ih€ H\n) :h>0,£ 

_SlaX{h>0} h€X m , 



T-\imF e {h) = F ao (h) 



oo otherwise, 



with x being the characteristic function. 

Proof. Consider an arbitrary sequence e„ — ¥ and h € X m . 
(i) (lim-inf condition) 

Let {h n } C X m such that h n — ^ h weakly in iJ 1 (SI) , then h n — > h strongly in L 2 (tt). 
Choose 5„ — > such that £ n /S n -> as n -> oo which immediately gives 



liminf / W e (hn) > liminf / W £ (h n ) = liminf / X {K > S n }. (5.3) 

™ Jn n J{h n >8 n } n Ja, 

Next we want to use J x{h > 0} < liminf J„ x{hn > 8 n }. Conversely assume 

liminf / X {K > S n } < [ X {h > 0}. (5.4) 
™ Jn J 

Then employing right-continuity of s t— > J n x{h > s } ( see 22|j Proposition 6.1) there 
must exist some 5, 5 > such that 

< liminf / X {h > 0} - x{K > $n} -S< liminf / X {h > 5} - x{K > S n }, 
71 Jn n Jn 

< liminf / x{h > 6 & h n < S n } < (t) liminf / \h- h n \ 2 = 0, 
11 Jn \oJ « Jn 
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where we used Chebyshev's inequality. This is a contradiction and thus by the previous 
assertions 

liminf / W e (h n )> [ X {h>0}. 
" Jn Jn 

(ii) (lim-sup condition) 

Define a recovery sequence by h n = a n h + e n where a n = (m— e n \fl\)/m. Then 
h n G X m and h n —> h even strongly in H 1 ^) and the following estimate holds 

limsup / W sn {K) = limsup ( / W ( 1 + —h] + [ W\l + —h] ) , 

n Jn n \J{h>0} \ £ n J J {h=0} \ £ n J ' 

< limsup / X {h > 0} + / W (1) = / x{h > 0} , 
n Jn J{h=o} Jn 

where we used that W(s) < 1 for s > 1 and W-^l) = 0. □ 

To prove the T-convergence to the sharp-interface model we can exploit the prop- 
erty that the behavior of gradient terms can be easily controlled. 

Theorem 5.3. For the family of energies (|5.1[) the T-limit is 

E 00 (h 1 ,h 2 )= [ ||V/ ll | 2 + i|V/ l2 | 2 + |$(l)|x{^2>M 
j n 

Proof. The gradient terms in E e are weakly lower semicontinuous with respect 
to weak convergence in iJ 1 (0) 2 . Together with Proposition 15.21 this gives the lim- 
inf inequality. On the other hand the gradient terms are continuous with respect to 
strong convergence in H 1 ^) 2 . Choosing the recovery sequence as in the proof of 
Proposition 15.21 one gets the desired lim-sup inequality. □ 

Now we want to deduce necessary conditions for minimizers of E^. We are 
especially interested in conditions at the points where the two-phase domain meets 
the one-phase domain. One problem is that one cannot immediately compute the 
Euler-Lagrange equations (L 2 -gradient) of the sharp- interface energy functional Eoo . 
This is due to J x{^2 > ^1} being only lower semicontinuous in the strong H 1 ^) 
topology, but not continuous nor diffcrcntiable. In fact directional derivatives of this 
part of the energy will almost surely be zero or infinite. Therefore we compute the 
directional derivative of E^ in another topology as follows: For ease of notation 
introduce 

uj = {x G ft : h 2 > hi} (5.5) 

and restrictions of h\ and hi to uj and f2 \ uj are called 

h\ := hi\u, hi := folw, h := /ii|o\w = MtA^, 

with boundary condition hi = hi = h on duj. We will now vary hi, hi and h but also 
uj. The formal calculation is restricted to smooth hi,h2,h and uj. Using these we can 
rewrite the energy using the restrictions as 



E^Qii.hi.h,^) = 
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where we define a' := 1 + a. A perturbation of r i— > (hi(r), h 2 (r), h(r), w(r)) can 
be parametrised using a diffcomorphism tp(o,r) : il —> SI with the property that 
lo(t) — {iI>(x,t) : x G w(0)}. In the same spirit define Jii(xq,t) — hi(ip(xo, t), t) as 
the pullback of hi by ip and similarly h 2 and h. The boundary conditions d T hi = 
d T h 2 — d T h on 8lo(t) translate into 



hi+tp- Vhi = h 2 +tp- V/i 2 = h + ip -Vh=: £. 



(5.7) 



Here we use the notation ip := d T ^p, hi := d r hi, h 2 := 5 r /i 2 and ft. := <9 T /i. Then 
using Reynolds transport theorem we get 



—Eoo(hi,h2,h,uj) 



crVftiV/li + V/l 2 V/l 2 
1 



a'VhVh 



di>J 



|V/n| 2 



|Vft 2 | 2 



a 1 



|V/i| 2 + 1*(1)| 



(n • VO 



Applying integration-by-parts and boundary conditions (|5 . T[) yields in one and two 
spatial dimensions, i.e. u C R and ui C R 2 the directional derivative 

( E °° +Mf hi+Xif (h 2 - hi)) , 
= - / [trA/i! + A 2 - Xi]hi + [Ah 2 +Ai]/i 2 - / [a'Ah + X 2 ]h 

J uj JQ,\lj 

f (VM 2 - \(Vh 2 f + y(Vft) 2 + |*(1)| + r, 2 (l + o a')] (n ■ i>) 

+ / [a(n-V/ii) + (n-V^-CT'(n-V/i)](£+(*-V>)). (5.8) 

where rj 2 := (t- V/ii) 2 = (t- V/i 2 ) 2 = (t-Wh) 2 . The expressions inside square brackets 
have to vanish independently, since the perturbations (hi, h 2 , h, £, ip) are independent. 

Remark 5.4. Above we added Lagrange multipliers Ai,A 2 to take care of the 
mass conservation. In one dimension there is no tangential contribution and hence 
7] = 0, whereas in two dimensions the contribution with rj 2 vanishes due to definition 
a 1 = a + 1. However, if a' could be choosen independently of a, there would be an 
extra contribution in that case. 

5.2. Existence and uniqueness of solutions. In this part we consider the 
sharp interface energy derived by T-convergence and study its minimizers. The idea 
of the proof is to show that for a minimizer the support of 



h := h 2 — hi 



(5.9) 



is a ball contained in 57, on which the solutions can be computed explicitly. The mini- 
mization itself is performed with masses m = (mi, m 2 ) held fixed. Further extensions 
of our proof and properties of the solutions are discussed in the end of this section. 



Definition 5.5. Let A c R™ a Borel set of finite Lebesgue measure, then the 
symmetric rearrangement of the set A is defined by A* = B s (0) with s such that 
fi(A) = fi(A*). The symmetric decreasing rearrangement of the characteristic function 
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is (xa)* ~ XA* ■ Now let f : R" — > R a Borel measureable function vanishing at 
infinity, then define the symmetric- decreasing rearrangement of f by 



f*(x) = / X{f >s }(x)ds. 



u 



Theorem 5.6. (Minimizer of sharp interface energy) Let = Br(0) and X = 
{(huhz) e X m (Cl) : {hi - h 2 )\ dn = 0} and energy 

E(hx,h 2 ) := jT ||Vfci| 2 + i|V/i a | a +\<f(l)\ X{h2>hl} dx. (5.10) 

Then using C( x ) := oi{s 2 — \x\ 2 ) + minimizers of E with mass (7711,7712) are 

hi= ° ((x - ,t ) + h, hi(x) = h2-((x-x ), (5-H) 
o + 1 

wi£/i constant xq £ £7 anc! r,a,h S R. Prescribing the mass (7711,7112) /ixes 7* a?i(i /i, 
whereas a is fixed by the contact angle (Neumann triangle) 

cr(V/ii) 2 + (V/i 2 ) 2 = 2|$(1)|, at\x\=r. (5.12) 

.For /arge masses ?ti 2 (|5.12p is not required and we get r = R, xq — in (|5.1ip . 



Proof. Using similar ideas as in [21[, we proceed as follows: 

Symmetry: For given (h\,h 2 ) £ X let h = (h 2 — hi) € /fg(f2) as in (|5.9[) . non- 
negative let 

l|V/i|| 

Using Cauchy-Schwarz inequality in L 2 (f2) the following energy estimate holds 

E{h u h 2 ) > (a + l)\\Vh 2 \\ 2 + o-\\\7h\\ 2 - 2<r\\\7h 2 \\ \\Vh\\ + |$(1)|/j({/i > 0}), 
= ||V/i|| 2 (A 2 +( 7(l-A) 2 )+|$(l)| A 7({/i>0}). 

Minimizing with respect to A gives the lower bound 

E(h!,h 2 ) > ^—\\Vh\\ 2 + |*(1)|a*({A > 0}) 
o + 1 

which is attained only if A = a/(a + 1) and Vh 2 is a multiple of Vh. Now let h* 
be the symmetric-decreasing rearrangement of h, then by virtue of the Polya-Szcgo 
inequality 



h* dx = h dx, ||V/i*||i3( R ») < ||V/i||i3(a»). 

JR™ 

We have the freedom to translate h* as long its sup port is contained in f2. Equality 
holds only if h is already symmetric-decreasing [28[. Now assume that h is not sym- 
metric decreasing, or V/i 7^ 0-/(0- + l)Vh 2 . Then we can reduce the energy by defining 
h\ and hi, by 

h 2 (x) := h*(x — Xq) + h, h\(x) := h* 2 (x) — h*(x — Xo), 

(7+1 
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and have Wh^ = XV h* and /j,{h > 0} = fi{h* > 0} so that 

E(h l7 h 2 ) > -^—\\Vh*\\ 2 + Ml)H{h* > 0}) = E(h* l7 h*). 
o + 1 

Note that by definition {h* > 0} = B s (0) =: uj* with s such that n{uo*) = n(h > 0). 
To check ((x) = a(s 2 — |a;| 2 ) + is now analogous to [2l|. One has to solve the Euler- 
Lagrange equation for the first variation of E given in (|5 ,8|) using standard methods. □ 



Corollary 5.7. Let X be as before and the sharp interface energy 

E^MM := j ^|V/u| 2 + i|V/i 2 | 2 + |$(l)|)da ; + jf y|V/i| 2 d.T 

(5.13) 

as in (|5.6|) for a' > arbitrary. Then the minimizers of (|5. 10|) and ()5.13j) m X are 
ideniica/. 



Proof. Since we have ft = on doj the estimates of the previous proof are valid if 

2 



the domain of integration is restricted to lo. By construction we have || V/i|j 2 2(m aj - ) > 



l|VA*||i aC n\ w . ) = 0.D 

Remark 5.8. Using the abbreviation c = |$(1)|(<7 + l)/cr we can easily compute 
the parameters r and a from the previous theorem and get 

9m 2 \ 1/4 ^2c^ 1/4 
Sid = I — — , au 



8c 7 ' ~ la V9m 2 



2\ 1/6 / 2 \ 1/3 

■ 2 c y \2rri2 / 

in one and two spatial dimension respectively. The contact angles are then 



n-V/iJ = ±-^Hl n-V/i 2 | =tV2c~— — . (5.14) 

and which actually holds in any spatial dimension. We also have 



n ■ V(/i 2 - hi) = V2c = ^- ^^ ^ $(1) (5-15) 

which can be compared with the appropriate boundary condition in (|4.25[) from the 
matched asymptotic expansion. 

Discussion and Outlook. We considered stationary solutions of systems of cou- 
pled thin-film equations for two-layer liquid films. After proving existence of station- 
ary droplet solutions, we used matched asymptotic analysis to derive a corresponding 
sharp-interface model in the limit when the thickness of the ultra-thin film in the 
dewetted region e — > 0, which then yields the equilibrium Neumann angles. We point 
out that our asymptotic analysis requires the inclusion of logarithmic switch-back 
terms for the asymptotic droplet solution, which should in principle also be needed 
for the limiting case of droplet solutions on solid substrates. 
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We then proved the existence and uniqueness of the sharp-interface model using 
the variational structure of the equations allowing us to formulate the problem as a 
minimization problem, for which we can study the limit e — > via T-convergence. In 
one spatial dimension on an interval both sharp- interface models are equivalent. In 
particular the contact angle of h from the matched asymptotic analysis in (|4.25[) is 
the same as the one from the T-convergence in (|5.15p . Since the recovery of h\,hi 
from h in both cases works via f)4. 231 14.24[) or (|5.11j) , the second contact angle agrees 
as well. We note that dimensions d > 1 one has to prove that the shape of the domain 
\hi — hi > 0} is a ball of a certain radius. Using symmetric decreasing rearrangement 
this property, and thereby existence and uniqueness of minimizers, could be proved. 

We expect that, as for thin films on solid substrate, the techniques of matched 
asymptotic analysis can be extended to the dynamic time-dependent problem. In 
particular, the derivation and study of the time-dependent sharp-interface model will 
also support the understanding of the energetic structure of the system of the coupled 
thin film model and should still be valid in the time-dependent problem, e.g. in 
the gradient flow structure of a sharp-interface model. This will be important in 
the study of dewetting regimes, dewetting rates and the stability properties of the 
evolving interfaces, as it was done previously for the dewetting liquid films from solid 



substrates, see e.g. [3JJ, |32 



As pointed out in the beginning of our study, mathematical theory for two-layer 
liquid flows leaves still many open questions and problems to be addressed. The 
present work can only be considered as a first step. Moreover, even considering only 
stationary solutions, we note that the general picture is much richer compared to the 
situation on a solid substrate, with energy structures leading to phase-inverted or 
more complicated patterns, and is be subject of our ongoing research. 
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